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High energy photon emission rate from matter created in 
Pb + Pb collisions at CERN SPS energies is evaluated. The 
evolution of matter from the initial state up to freeze-out has 
been treated within the framework of (3+1) dimensional hy- 
drodynamic expansion. We observe that the photon spectra 
measured by the WA98 experiment are well reproduced with 
hard QCD photons and photons from a thermal source with 
initial temperature ~ 200 MeV. The effects of the spectral 
changes of hadrons with temperature on the photon emission 
rate and on the equation of state are studied. Photon yield 
for An + An collisions at RHIC energies is also estimated. 

PACS: 25.75.+r;12.40.Yx;21.65.+f;13.85.Qk 



Ultra-Relativistic collisions of heavy nuclei have 
brought us within reach of creating and studying vari- 
ous aspects of quark-gluon plasma (QGP), which so far 
was believed to exist in the microsecond old universe or 
possibly in the cores of neutron or quark stars. We are at 
a very interesting situation in this area of research where 
the Super Proton Synchrotron (SPS) era has drawn to a 
close and the first results from the Relativistic Heavy Ion 
ColHder (RHIC) have started to appear. Already from 
the results of the Pb run at the SPS quite a few of the 
signatures of QGP, e.g., J/^ suppression, strangeness 
enhancement etc., are reported to have "seen" unmistak- 
able hints of the existence of QGP . Electromagnetic 
probes, viz., photons and dileptons have long been recog- 
nized as the most direct probes of the collision ||^ . Owing 
to the nature of their interaction they undergo minimal 
scatterings and are by far the best markers of the entire 
space-time evolution of the collision. 

The single photon data, obtained from Pb-Pb colli- 
sions at CERN SPS reported by the WA98 Collabora- 
tion have been the focus of considerable interest in 
recent times While in Ref. the contribution 

from hard photons and the effects of transverse expansion 
on the thermal photons were neglected, in the ther- 
mal contribution was not taken into account. Again, in 
1^,0 the thermal shift on hadronic masses was neglected, 
whereas, in |^ both the the intrinsic kx distribution of 
partons and thermal shift of hadronic masses were ig- 
nored. In this letter we present the results of an analysis 
of the photon spectra in a realistic framework using a 
reasonable set of parameters and consistently taking into 
account all the effects mentioned above. We emphasize 
the effects of in-medium modifications of hadrons on the 
photon spectra considering the fact that as yet it has 
not been possible to explain the observed low-mass en- 



hancement of dileptons measured in the Pb+Au as well 
as S-|-Au collisions at the CERN SPS in a scenario which 
does not incorporate in-medium effects on the vector me- 
son mass (see ||^ for a review). 

Let us first identify the possible sources of "excess" 
photons above those coming from the decays of pseu- 
doscalar tt*^ and r] mesons, as provided by the data. 
Firstly, one has the prompt photons coming from the 
hard collisions of initial state partons in the colliding 
nuclei. These populate the high transverse momentum 
region and can be estimated by perturbative QCD. The 
thermal contribution depends on the space-time evolu- 
tion scenario that one considers. In the event of a de- 
confinement phase transition, one first has a thermalized 
QGP which expands and cools, reverts back to hadronic 
matter, again expands and cools and eventually freezes 
out into hadrons most of which are pions. Photon emis- 
sion in the QGP occurs mainly due to QCD annihilation 
and Compton processes between quarks and gluons. In 
order to estimate the emission from the hadronic matter 
we will consider a gas of light mesons viz. tt, p, oj and rj. 

It has been emphasized that the properties of vector 
mesons may change appreciably because of interactions 
among the hadrons at high temperatures and/or densi- 
ties (see the reviews |^-|l^). This modifies the rate of 
photon emission as well as the equation of state (EOS) 
of the evolving matter. Among various models for vector 
mesons available in the literature [p^, we examine two 
possibilities for the hadronic phase in this paper: (i) no 
medium modifications of hadrons, and (ii) the scenario 
of the universal scaling hypothesis of the vector meson 
masses [ pT| . In principle, we can think of a third sce- 
nario (iii) the large coUisional broadening of the vector 
mesons H . Both (ii) and (iii) can reproduce the enhance- 
ment of the low-mass dileptons measured by CERES Col- 
laboration at CERN SPS, but the scenario (iii) has been 
found to have a negligible effect on the emission rate of 
photons |l3|. The effect of temperature dependent mass 
as described in case (ii) has also been incorporated in the 
EOS of the hadronic matter undergoing a (3-1-1) dimen- 
sional expansion. We will see that the resulting photon 
spectra reproduce the experimental data quite well. 

There is still substantial debate on the order of the 
phase transition as well as the value of the critical tem- 
perature (Tc). To address this aspect we will also con- 
sider a scenario where the system begins to evolve from 
a high temperature phase where all the hadronic masses 
approach zero (pion mass is fixed at its vacuum value). 
As the system expands and cools, the hadrons acquire 
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masses (as in case (ii) above) till freeze out. Incorpora- 
tion of medium modified masses and the EOS in this case 
also provides a reasonable explanation of the data. 

We start with the direct QCD photons originating from 
hard scattering of partons embedded in the nuclcons of 
the colliding nuclei in the very early stages of the colli- 
sion. We will see later that they make a significant con- 
tribution to the total photon yield. These are estimated 
using perturbative QCD as 
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where TAA{b) is the nuclear thickness at impact param- 
eter b. Its value at = is taken as 220/fm^ g]. app 
includes the pp cross-section for Compton and annihi- 
lation processes among the partons. At SPS energies 
one should include the effects of intrinsic kx distribution 
of partons to account for the fact that the collid- 
ing partons might have some initial transverse momenta 
with respect to the incoming hadrons (analogous to the 
Fermi motion of nuclcons in the nucleus). This leads to 
substantial enhancement in the photon spectra [^. In 
practice such an effect is implemented by multiplying 
each of the parton distribution functions appearing in 
the right hand side of the above equation by a Gaussian 
function of the type /(fc^) = exp[— fc^/(fc^)]/7r(fc^) and 
integrating over (Pkr- We use CTEQ5M partons and 
(k^) = 0.9 GeV2 g for evaluating the hard QCD pho- 
tons (•\/< k^ > = 0.8 GeV is taken in |6|). The photons 
from hard QCD processes have been used to normalize 
the p-p data. The higher order effects has been taken into 
account through a K-factor ~ 2. An enhanced produc- 
tion in A-A collisions compared to p-p will presumably 
mark the presence of a thermal source . 

The lowest order processes contributing to hard ther- 
mal photon emission from quark gluon plasma are the 
QCD Compton and annihilation processes It has 

been shown recently that the two-loop contribution 
leading to bremsstrahlung and qq annihilation with scat- 
tering is of the same order as the lowest order processes.. 
The total rate of emission per unit four- volume at tem- 
perature T is given by 
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where Jt — 4.45 and Jl 
'g' is given by, 

5^/47r = as ^ 



-4.26. The QCD coupling, 
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for two quark flavours ||l8| . For photon energies in the 
range 1 < E (GeV) < 5, the static emission rate given 
by Eq. (|) at T ^ 200 MeV is about an order of mag- 
nitude larger than the rate due to Compton and anni- 
hilation processes |10]. We will see below that within 



the present framework the space time integrated photon 
yield from quark matter is less than that from hadronic 
matter due to the smaller life time of the QGP phase as 
a result of a moderate value of the initial temperature 
considered. Therefore, the overall thermal photon yield 
remains largely unaffected by the two loop contribution. 

To estimate the photon yield from the hadronic matter 
(HM) (see first of jl^l), we have considered the reactions, 
7rp —> 7r7, TTTT —>■ p7, TTvr 7r?7 —> tt ^ and the 

decays p nn j and lo tt"/. The invariant ampli- 
tudes for all these processes are given in Refs. [|l9|. In the 
present work we have also considered photon production 
due to the process tt p ^ ai — > tt 7 ||l^ . 

In Ref. ||l^ we have studied the effects of spectral 
changes of hadrons on the electromagnetic probes in de- 
tail. It was observed that the gauged linear and non- 
linear sigma models and the model with hidden local 
symmetry do not show any appreciable effect on photon 
emissions. In the Walecka model, the universal scaling 
hypothesis for the vector meson masses as well as the 
large collisional broadening of vector mesons produce a 
large enhancement in low mass dileptons. However, the 
photon emission rate does not suffer substantial medium 
effect in the latter case, since the spectral function is 
smeared out. Nevertheless, the scaling hypothesis with 
particular exponent A = 1/2 (called the Nambu scaling 
in Jll| ) has been seen to enhance photon emission among 
the others. 

To consider the effect of the spectral modifications of 
hadrons we adopt two extreme cases: (i) no medium 
modifications of hadrons, and (ii) the scaling hypothe- 
sis with A = 1/2. In case (ii), the parametrization of 
in- medium quantities (denoted by *) at finite T is 
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where V stands for vector mesons, fv is the coupling 
between the electromagnetic current and the vector me- 
son field and luq is the continuum threshold. Mass of the 
nucleon also varies with temperature as Eq. (^). Note 
that there is no definite reason to believe that all the in- 
medium dynamical quantities are dictated by a single ex- 
ponent A. This is the simplest possible ansatz (see [|l^ for 
a discussion) . The effective mass of ai is estimated by us- 
ing Weinberg's sum rules . We have seen earlier that 
the baryon chemical potential has a small effect on the 
photon yield [Q. Moreover, in the central rapidity region 
the entropy per baryon is quite large '--^ 40 — 50 pl| , p2| . 
Thus the finite baryon density effects are neglected here. 

We will assume that the produced matter reaches a 
state of thermodynamic equilibrium after a proper time 
~ 1 fm/c If a deconfined matter is produced, it 

evolves in space and time till freeze-out undergoing a 
phase transition to hadronic matter in the process. The 
(3-1-1) dimensional hydrodynamic equations have been 
solved numerically by the relativistic version of the flux 
corrected transport algorithm [e3, assuming boost in- 
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variance in the longitudinal direction |2^] and cylindrical 
symmetry in the transverse plane. The initial temper- 
ature Ti can be related to the multiplicity of the event 
dN/dy by virtue of the isentropic expansion as p5|] , 
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where Ra is the initial radius of the system, r.^ is the 
initial thermalization time and Uk — (tt^/OO) gk, Qk be- 
ing the effective degeneracy for the phase k (QGP or 
hadronic matter). The bag model EOS is used for the 
QGP phase. gniT), the statistical degeneracy of the 
hadronic phase, composed of tt, p, uj, 77, ai and nucleons is 
a temperature dependent quantity in this case and plays 
an important role in the EOS As a consequence the 
square of sound velocity, c^^ SH)idgH /dT) + 3] < 

1/3, for the hadronic phase, indicating non- vanishing in- 
teractions among the constituents (see also [§6|). The 
hydrodynamic equations have been solved with initial en- 
ergy density, e{Ti,r) [Q, obtained from Ti through the 
EOS. We use the following relation for the initial velocity 
profile which has been successfully used to study trans- 
verse momentum spectra of hadrons [^,^ and also 
photons 
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For our numerical calculations we choose ^ = 1 and sen- 
sitivity of the results on vq will be shown. It is observed 
that the results do not change substantially with reason- 
able variation of the parameter S for a given value of vq. 
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FIG. 1. Transverse momentum distribution of thermal 
photons with and without medium effects. 

For central collisions of Pb nuclei at 158 AGeV at 
the CERN-SPS, we assume that QGP is produced at 
Ti—1 fm/c which expands in (3-1-1) dimension and under- 
goes a first order phase transition to hadronic matter at 
rc=160 MeV. Taking diV/d?/=700 (last of Ref. §), and 
gk = gQGP=-i7 for a two-flavour QGP, the initial tem- 
perature Ti comes out as 196 MeV. In a first order phase 



transition one has a mixed phase of coexisting QGP and 
hadronic matter which persists till the phase transition 
is over. Thereafter the hadronic matter expands, cools 
and freezes out at a temperature, Tf and radial velocity, 
vf. The sum total of the photon yields from the QGP 
phase, the mixed phase and the hadronic phase consti- 
tutes the thermal yield. The values of (T/,w/) should 
in principle be obtained from the analysis of hadronic 
spectra. However, at present it is not possible to do so 
without ambiguity. In Ref. it is shown that the ex- 
perimental data from Pb + Pb collisions allow values of 
{Tf,vl) ranging from (180 MeV,0) to (120 MeV,0.7). In 
the present work we take Tf — 120 MeV and the initial 
velocity profile of Eq.(H) with vq as a parameter. 
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FIG. 2. Total photon yield in Pb -I- Pb collisions at 158 
A GeV at CERN-SPS. The theoretical calculations contain 
hard QCD and thermal photons. The system is formed in the 
QGP phase with initial temperature Ti — 196 MeV. 

In Fig. ^ we show only the thermal photon spectra 
originating from quark matter (QM = QGP -I- QGP part 
of mixed phase) and hadronic matter (HM = hadronic 
part of mixed phase + hadronic phase) with vq = 0. The 
solid and dashed lines correspond to the case (i) and (ii) 
respectively. The case with large collisional broadening 
shows no deviation from (i). At low px the difference 
between (i) and (ii) is negligible because in this region of 
phase space most of the photons are emitted from the late 
stage of the evolution where the in-medium effects are 
small. The increased photon yield at large px is caused 
by the enhancement in the Boltzmann factor due to the 
reduction in meson (particularly, p) masses. However, in 
the total photon emission this difference of thermal pho- 
tons at high pt is masked by the hard photon contribu- 
tion. The thermal photon yield with hadronic mass vari- 
ation due to Walecka model or Brown-Rho scaling |l|] 
(A = 1/6 in Eq. (^) will lie in-between the two curves in 
Fig. 0. 

In Fig. 1^, results for the total photon emission is shown 
for three different values of the initial transverse velocity 
with medium effects as in case (ii). All the three curves 
represent the sum of the thermal and the prompt photon 
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contribution which includes possible finite fc^ efi'ects of 
the parton distributions. The later, shown separately by 
the dot-dashed line also explains the scaled pp data from 
E704 experiment ||2^ . We observe that the photon spec- 
tra for the initial velocity profile given by Eq. (|) with 
vq — 0.2 explains the WA98 data reasonably well. A sim- 
ilar value of the initial radial velocity has been obtained 
in iQ from the analysis of transverse momentum spectra 
of hadrons and photons (see also Q). It is found that a 
substantial fraction of the photons come from mixed and 
hadronic phase. The contribution from the QGP phase 
is small because of the small life time of the QGP (~ 1 
fm/c). 
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FIG. 3. Total photon yield in Pb -I- Pb collisions at 158 
A GeV at CERN SPS. The theoretical calculations contain 
hard QCD and thermal photons. The system is formed in the 
hadronic phase with the hadronic masses approaching zero at 
initial temperature Ti — 205 MeV. 

The last statement together with the current uncer- 
tainty of the critical temperature Tc poses the fol- 
lowing question: Is the existence of the QGP phase es- 
sential to reproduce the WA98 data? To study this 
problem, we have considered two possibilities: (a) pure 
hadronic model without medium-modifications, and (b) 
pure hadronic model with scaling hypothesis according 
to Eq.(^. In the former case, Tj is found to be ~ 250 
MeV for t; = 1 fm/c and dN/dy — 700, which appears 
to be too high for the hadrons to survive. Therefore this 
possibility should be excluded. On the other hand, the 
second case with an assumption of Ti = Tc (which is just 
for simplicity) leads to Ti ~ 205 MeV, at = 1 fm/c, 
which is not unrealistic. In this case, the hadronic system 
expands and cools and ultimately freezes out at T/=120 
MeV. The masses of the vector mesons increase with re- 
duction in temperature (due to expansion) according to 
Eq.(^. The results of this scenario for three values of 
the initial radial velocity including the prompt photon 
contribution are shown in Fig. |^. The experimental data 
are well reproduced for vanishing initial transverse veloc- 
ity also. This indicates that a simple hadronic model is 
inadequate. Either substantial medium modifications of 



hadrons or the formation of QGP in the initial stages is 
necessary to reproduce the data. It is rather difficult to 
distinguish between the two at present. 

We now show our prediction of the photon yield in 
central collisions of Au nuclei at 200 A GeV at RHIC 
in Fig. ^ for scenario (ii) with vq = 0. Enhancement 
due to intrinsic motion of the partons has been ignored 
at RHIC, because at higher beam energies such effects 
are small The initial temperature and thermaliza- 
tion time are taken as 300 MeV and 0.5 fm/c respec- 
tively iQ . It is seen that up to pt = 2 GeV, most of the 
photons are emitted from the thermal source. We also 
observe that photons from quark matter make a sub- 
stantial contribution to the thermal yield (dash-dotted 
line in Fig. ^), because of the higher initial temperature 
and larger lifetime of the QGP phase realized at RHIC 
compared to SPS. The effect of the transverse expansion 
on the QGP phase is small. However, photons from the 
hadronic phase (short dashed line in Fig. ||), particularly 
during the late stage of the evolution, receive a large kick 
due to the radial expansion and consequently populate 
the high px region, as shown in Fig. 0. 



10 

lo" 

10"' 
10"' 

;io"' 

10"'' 
10"' 
10"' 
10"' 



QM 
HM 

Thermal (QM+HM) 
Total (Tbermal+Hard) 




RHIC 

Au+Au @ 200 A GeV 
T|=300 MeV 
Ti=0.5 fm/c 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

p^(GeV) 

FIG. 4. Photon spectra at RHIC for Au + Au collisions. 

In summary, we have evaluated the high energy pho- 
ton yield in Pb -|- Pb collisions at CERN SPS energies 
with two different initial conditions. In the first sce- 
nario, we start with the assumption of the formation of 
a QGP phase at Ti ~ 196 MeV and then the system con- 
tinues through mixed phase and hadronic phase before 
freeze-out. In the second scenario, we assume a chirally 
symmetric phase where the hadronic masses approach 
zero at a temperature ~ 205 MeV and then the system 
evolves towards freeze-out. The effects of the variation of 
hadronic masses on the photon yield have been consid- 
ered both in the cross section as well as in the EOS. The 
photon spectra reported by the WA98 collaboration are 
well reproduced in both cases. We thus conclude that the 
thermal photon spectra resulting from the Pb -I- Pb col- 
lisions at CERN SPS energies are emitted from a source 
with an initial temperature ~ 200 MeV. A similar value 
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of Ti is also obtained from the analysis of photon and 
dilepton data from CERN SPS At RHIC energies 

the total yield of photons increases by an order of mag- 
nitude compared to SPS with a substantial contribution 
from QM. 

In spite of the above encouraging situation, a firm con- 
clusion about the formation of the QGP at SPS neces- 
sitates a closer look at some pertinent but unsettled is- 
sues. In the evaluation of the hard photon contribution 
the value of the K-factor and the intrinsic fcy distribu- 
tion of partons are adjusted so as to reproduce the scaled 
p-p data of E704 collaboration. However, it is extremely 
important to know quantitatively the contribution from 
the hard processes. The possibility of A and kx depen- 
dence of the K-factor |^^, the energy loss of fast partons, 
the shadowing of the structure functions and the broad- 
ening of transverse momentum distribution of partons 
are some of the vital issues related to photon produc- 
tion in nucleus-nucleus collisions which have not been 
considered in the present work. Again, the assumption 
of complete thermodynamic equilibrium for quarks and 
gluons may not be entirely realistic for SPS energies; lack 
of chemical equilibrium will further reduce the thermal 
yield from QGP. We have assumed = 1 fm/c at SPS 
energies, which may be considered as the lower limit of 
this quantity, because the transit time (the time taken 
by the nuclei to pass through each other in the CM sys- 
tem) is ~ 1 fm/c at SPS energies and the thermal system 
is assumed to be formed after this time has elapsed. In 
the present work, when QGP initial state is considered, 
we have assumed a first order phase transition with bag 
model EOS for the QGP for its simplicity, although it is 
not in complete agreement with the lattice QCD simula- 
tions [^o| . However, it is difficult to distinguish among 
different EOS with the current resolution of the photon 
data. As mentioned before, there are uncertainties in the 
value of Tc , a value of Tc ~ 200 MeV may be consid- 
ered as an upper limit. Moreover, the photon emission 
rate from QGP given by Eq. evaluated in Refs. [ p^Jl?) 
by resumming the hard thermal loops is strictly valid for 
g « \ whereas the value of g obtained from Eq. is 
~ 2 at T ~ 200 MeV. At present it is not clear whether 
the rate in Eq. (Q) is valid for such a large vaule of g or 
not. 
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